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We extend a fictitious-cavity approach to calculate the Casimir eflect for cavities bounded by 
flat anisotropic materials. We calculate the energy, force and torque in terms only of the optical 
coefficients of the walls of the cavity. We calculate the Casimir effect at zero and finite temperature 
for some simple systems. As a non trivial application, we calculate the torque between a semi-infinite 
anisotropic plate and an anisotropic film. We study the effect of the film thickness in the torque 
and find an optimal width that maximizes the torque. 
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I. INTRODUCTION 

An electromagnetic mode of frequency w within an electromagnetic cavity is analogous to a harmonic oscillator 
with a quantized energy spectrum given by semi-integer multiples n -f- 1/2 of the energy quantum hui, where the 
integer n is the occupation number of the mode, i.e., the number of photons in the corresponding state, Tiu is 
the quantized energy of each photon and Tiuj/2 is the ground state energy arising from the quantum nature of the 
electromagnetic field and its zero point fluctuations. As the frequency of the electromagnetic modes depends on the 
geometry of the cavity, the zero point fluctuations may not be simply disregarded as a constant contribution to the 
energy and should be accounted for in calculations of the total energy of the system. A simple vacuum cavity may 
be produced by positioning two flat conducting plates parallel to each other a small distance L apart, leading to 
a first quantization of the electromagnetic field within the cavity. In 1948, Casimir [T] predicted that the vacuum 
electromagnetic energy due to the quantum fluctuations of the field would depend on L and therefore a force, which 
turns out to be attractive, would act on each plate. This Casimir force also may be conceived as originated from the 
difference between the radiation pressure due to the fluctuating electromagnetic modes outside the plates and the 
modes within [5]. Irrespective of its interpretation, the Casimir force has its origin in the linear momentum carried by 
the radiation field. Nevertheless, the radiation field carries angular momentum beyond linear momentum. Therefore, 
if the plates are optically anisotropic, a Casimir torque might develop besides the Casimir force. In fact, the transfer 
of angular momentum of polarized light to a macroscopic birefringent medium, resulting in a torque, has been known 
for a long time ^ . The Casimir torque may be interpreted as arising from the orientational dependence of the vacuum 
electromagnetic energy, specifically from the dependence on the relative orientation of the optical axes of the plates. 
The resultant torque tends to align the optical axes along the configuration that minimizes the vacuum energy. 

Recently, the Casimir effect has received considerable attention for its possible technological applications, besides 
the fact that experimental studies have attained the necessary accuracy to test in detail the theoretical predictions 
[4HT7]. Therefore, theories about the Casimir effect that account realistically for the properties of actual materials 
have become indispensable. The study of vacuum forces between real materials was pioneered by Lifshitz [18] . who 
considered two semi-infinite homogeneous and isotropic non-spatially dispersive dielectric slabs, whose fluctuating 
currents were the sources of the fluctuating electromagnetic field and whose correlations were related to the dielectric 
response of the materials. 

In 1972, Parsegian and Weiss derived an expression in the non-retarded limit for the interaction energy between 
two semi-infinite dissipationless dielectric anisotropic materials [19' following a method of surface mode summation. 
Barash derived an expression which included retardation and dissipation effects [20 employing an auxiliary system 
[2T] first introduced for isotropic systems. The solutions of Maxwell's equations for the field in inhomogeneous and 
absorbing media were expanded in terms of the orthogonal solutions of Maxwell's equations for a non-dissipative 
auxiliary system in which the frequency dependence of the dielectric function is only parametric. 
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The use of an auxiliary system was further developed in physical terms by Kupiszewska [22] in a calculation of 
Casimir forces for lossy and dispersive isotropic dielectrics in the case of one dimensional propagation of electromag- 
netic waves. The problem of quantizing a dissipative system is attacked by accounting both for the dynamics of the 
vacuum modes and of the atomic dipoles to which they couple and which make up the material, together with a 
thermal reservoir in which the atomic radiators dissipate the absorbed energy. That formalism was extended by van 
Enk 23 to obtain the torque between anisotropic materials in the ID case. In his work, van Enk calculated the 
torque starting from the flux of the spin angular momentum of the electromagnetic field. 

Numerical calculations have also been performed for materials with a small anisotropy using Barash's results and it 
has been shown that the torque may be large enough to be experimentally measurable in several novel experimental 
configurations j24[. By a similar technique, the Casimir energy between anisotropic dielectric plates |25j . and between 
a plate with anisotropic magnetic response and another with anisotropic dielectric response |26| have been calculated 
and analytical approximate expressions for the torque and force were obtained in the retarded limit. Kenneth and 
Nussinov have also calculated the Casimir energy between parallel plates made up of arrays of wires aligned along 
different directions [27] using a path integral technique [28l |29] . 

In the works described above, specific models of the dielectric properties of the plates were assumed from the onset 
in order to derive expressions for the Casimir force and torque. However, recent works |30H36) have shown that if 
the theory is set up in terms of the refiection coefficients of the media, or equivalently, in terms of their exact surface 
impedance [37l|38], it is possible to decouple the calculation of the Casimir force from the calculation of the dielectric 
response of the materials. 

Lambrecht et al. have extended their scattering approach [301121] to corrugated systems [33] ■ Moreover, they have 
argued that the resulting formula for the Casimir energy has a wider range of applicability and may be used to study 
other anisotropic mirrors. Their formula has been used to evaluate numerically the effects of corrugation on the 
Casimir force [ID] and torque [H] and to calculate the Casimir force between anisotropic metamaterials [32]. 

Mochan et al. |32H36) have argued that in thermal equilibrium, all of the properties of the radiation field within a 
cavity are completely determined by the optical reflection amplitudes of the walls. Indeed, whenever a photon reaches 
the surface of a wall of the cavity it may be coherently reflected with an amplitude described by the optical coefflcients 
of the wall. Otherwise, it would be transmitted into the wall to be either absorbed, exciting the material degrees of 
freedom of its constituents, or transmitted across the wall and into the surrounding vacuum to be lost forever. The 
probability of these processes is again determined by the optical coefficients of the wall, and given by their squared 
modulus. In thermodynamic equilibrium, detailed balance implies that whenever a photon is lost, an equivalent 
photon is incoherently injected back to the cavity, so that both, the coherently and incoherently reflected photons are 
determined by the reflection amplitudes alone. Therefore, the radiation fleld within a real cavity would be identical 
to the field within any cavity that has walls with the same optical properties. This fact allowed the construction of 
fictitious disipationless systems which can be simply treated quantum mechanically to obtain the vacuum energy and 
force for cavities with arbitrary walls. Thus, expressions for the Casimir force obtained from the electromagnetic stress 
tensor can be applied to semi-infinite or finite, homogeneous or layered, local or spatially dispersive, transparent or 
opaque systems through a simple substitution of the appropriate optical coefficients. This formalism has allowed the 
calculation of the Casimir force between photonic structures |43j . non-local excitonic semiconductors ,44) . non-local- 
plasmon-supporting metals with sharp boundaries [32H361 H5] . and between realistic spatially dispersive metals with 
a smooth self-consistent electronic density proflle [46]. With a few modifications, it has also been employed for the 
calculation of other macroscopic forces, such as those due to electronic tunneling across an insulating gap separating 
two conductors [17] . 

The relative simplicity of the formalism developed in |32H36j has allowed its generalization to anisotropic systems 
[48j . In Ref. |48j a new derivation of the Casimir torque within ID cavities with walls made up of arbitrary materials 
characterized only by their anisotropic optical coefficients was presented. By ID cavity we mean one in which the 
field is constrained to propagate only along one direction, namely, the normal to the surface of the cavity walls. In 
the present paper we generalize this formalism to 3D cavities with anisotropic walls. We calculate the Casimir force 
from the electromagnetic stress tensor. A simple integration over the separation distance L yields then the vacuum 
energy. The torque is then calculated by taking the derivative respect to the angle 7 between the optical axes of 
the plates. As our formalism is based on the calculation of the force, which is a directly observable quantity, it 
avoids the cumbersome singularities that plague other approaches. Furthermore, our results are written directly in 
terms of the optical coefficients of the walls of the cavity about which we make no assumption. Thus, they can be 
applied immediately to manifold systems such as insulating or conducting anisotropic slabs either dissipationless or 
dissipative, to semiinflnite or flnite walls and to homogeneous or structured materials. We test the validity of our 
approach by reproducing some known results [23, 25, 27, 30 36, 42 and we show its versatility by applying it to some 
previously unexplored systems. 

The structure of the paper is the following: In section [IT] we develop our formalism in order to arrive at expressions 
for the Casimir energy between arbitrary anisotropic plates. The use of an effective cavity allows our results to be 
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FIG. 1: (a) Vacuum cavity V of width L bounded by two arbitrary anisotropic material slabs (1 and 2) with surfaces at z\ 
and Z2 and anisotropic reflection amplitudes r'^" {a — 1,2, fijV — s,p). (b)Fictitious system made up three empty regions 
/, //, and III, bounded by perfect mirrors at zo and 23 and with infinitely thin sheets at zi and Z2 with identical reflection 
amplitudes r^'' to those of the real system and corresponding transmission amplitudes t'^'^ . 

applicable to arbitrary slabs at any temperature T. In section[TTT]we specialize our results to semi-infinite local uniaxial 
media and we apply our results to a calculation of the energy and torque for an idealized uniaxial system consisting of 
anisotropic mirrors that are perfectly conducting along one direction and perfectly insulating along the perpendicular 
directions. These calculations are performed for temperatures T = and T ^ 0. We also consider systems with a 
finite frequency-dependent conductivity along the optical axis. As a further application of our formalism, in section 
|IV| we calculate the torque of a system consisting of two conducting plates with an anisotropic effective mass tensor, 
where one of the plates is a film with finite thickness d, and we find there is an optimum value of d that maximizes 
the torque. We also calculated the same system considered by Munday et al. [24] consisting of a semi-infinite slab of 
BaTi03 and a thin film of calcite. Finally, in section [V] we present our conclusions. 

II. THE EFFECTIVE CAVITY APPROACH 

Consider the setup shown in Fig. [l|a). The slabs represents arbitrary media. According to H51 149 ) . in 

thermodynamic equilibrium the properties of the radiation field within the cavity V are completely determined by the 
geometry of the cavity, characterized by L, and by the 2x2 reflection amplitude matrices r^'^ of each slab (a = 1,2) 
coupling z^-polarized incident light to /i-polarized reflected light {fj,, v = s,p). Any relevant property of the material is 
completely accounted for through its optical coefficients. Thus, the electromagnetic radiation within the real cavity 
V must be identical to that within a fictitious cavity V = // bounded by infinitely thin sheets at zi and Z2, provided 
their reflection amplitudes r^'^ are chosen to match those of the walls of the real cavity V. The transmission amplitudes 
f^" of the inflnitely thin sheets are conveniently chosen in order to guarantee energy conservation with no absorption 
whatsoever of electromagnetic energy. Thus, there is no absorption in the flctitious system, there is no excitation of 
material degrees of freedom and the normal modes of its electromagnetic fleld form a complete orthogonal basis of the 
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corresponding Hilbert space. As a consequence, one is allowed to use well developed quantum-mechanical procedures 
for the calculation of the field properties without the requirement of a microscopic model of the material. 

The field modes may be quantized and counted by choosing suitable boundary conditions. For example, we can 
can add perfect mirrors far away from the walls of the real cavity (Fig. [ijb)), at zq and 2:3. These quantizing mirrors 
produce a field that mimics the incoherent radiation back into the cavity that is responsible for maintaining a detailed 
balance and thus the thermodynamic equilibrium. 

Consider now a single wave of frequency lu with wave-vector projection Q along the interface. Without loss of 
generality, we choose x — 2; as the plane of incidence, so the electric and magnetic fields are given by 



and 
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where q — lu/c is the free-space wavenumber, (t>''{z) and 4>''{z) are the p and s polarized components of a spinorial 
normalized wave-function 
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where C^^ are constant coefficients within each region A = /, //, /// corresponding to /x-polarized light moving 
towards the right (C — r) and left {( — I), and k = \J\ijp' jc? — Q^) is the wave-vector component perpendicular to 
the interface. 

We integrate the energy density. 
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to obtain the total electromagnetic energy 
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in the limit Lj,Ljjj -> cjo, where A is the area of the plates and ||Ca|P = '^^c I^A^I^- ^''^ ^^'^ same limit the 
normalization condition imposed on the wave-function simplifies to 



l^{Lj\\Ci\\^ + Lui\\Ciii\\^). 
Notice that most of the energy lies in the large fictitious regions I and III, so we may identify 

A\£o\' 



U = 
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and solve for the amplitude |£oP = ?>t^U/A. To obtain the force on the slab 2, we calculate the stress tensor 
Tij ^ (l/87r)Re[£:i£;* + B^B* - {\E\'^ + \B\'^)5ij/2] at an arbitrary position z within the cavity, 
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By applying boundary conditions at zq and 23, we obtain for a given value of Q a discrete set of mode frequencies 
LUn and corresponding perpendicular components fc„ of the wave- vector. Each of these modes contributes to the stress 
tensor a quantity similar to that in Eq. ([8|), so that 
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where we have substituted the energy Un in terms of the equilibrium occupation number — f{ujn) — coth(/3?icj„/2)/2 
of a photon state with quantized energy hujn at temperature ksT = with ks the Boltzmann's constant. The 
sum over states may be rewritten in terms of the tensorial Green's function G(z, z') with components 
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for the ID Helmholtz equation 



(11) 



Here, k = k + irj, (77 > 0), with the understanding that the hmit — 0+ is to be taken at the end of the calculation. 
Using the identity Im(^^ — /c^)"^ = —T:6{k^ — fc^), we can replace the sum Kn by the integral 
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where / — f{uj) and 
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(13) 



plays the role of a local density of states at z (number of states per unit length and per unit fc^). 

The solution of ( 11 ), subject to the appropriate boundary conditions, may be written in terms of the solutions u(z) 
and v(z) of the ID Helmholtz equation that satisfy the boundary conditions on the right and left side of the system, 
respectively. 



G{z,z') = u(z)[u'(z')-v'(z')v-i(z')u(z')] ^e{z-z') 

-v(z) [v'(z') - u'{z')u-\z')v{z')] '^0{z' - z). 



(14) 



where 9 denotes the Heaviside unit step function. Here, u(z) and v(z) are 2x2 matrices with matrix elements u^{z) 
and v'^^z), A = 1,2 denotes the two independent spinorial solutions of Helmholtz equation, /i — s,p denotes their 
correspondent s and p components and u'(z) and v'(z) denote the derivatives of u(z) and v(z) with respect to their 
argument. We remark that we first introduced a similar expression for the spinorial Green's function in Ref. |48j . 
where it was used to obtain the angular momentum flux within a ID cavity. 

The solutions u(z) and v(z) may be written within the cavity in terms of the reflection coefficients of the plates 
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(a = 1,2). 



These are defined through 
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where we define the unnormalized spinors ^^id — s,p) as the amplitudes of the incident and reflected fields 

at the surface of plate 2 through 



and 
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The other components of the electromagnetic field may be obtained from Eqs. (17 1, (18 1 through Maxwell's curl 
equations. The matrix elements r^'^ of plate 1 are similarly defined. Thus, we may write. 



and 



where I is the 2x2 unit matrix. 
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Substitution of (19) and (20) in (14), (13) and (12) yields 
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where Va = det Va and 



A^det (I- e^^^'-^rirs). 

Note that A = yields the dispersion relation of the lossy modes of the real cavity. 
Smnining equation (21) over Q we finally obtain 



(22) 



he 

2^ 



(fQ Re 



(23) 



where we assumed Born-von Karman periodic boundary conditions along the surface of area A — oo to replace 
EQ...^A/(27r)2/d2Q... 

The flux of linear momentum —T^z in the fictitious cavity is the same as in the real cavity between slabs 1 and 2. 
To obtain the force on slab 2, we have to subtract the fiux in the real system between the slab and infinity. This can 
can be obtained following the same derivation given above, but replacing the slab 1 by the complete system made up 
of slabs 1, the cavity V and slab 2, and replacing slab 2 by empty space. The result is identical to equation (23), but 
substituting r2 — > 0. Thus, the total force per unit area of slab 2 is 
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The derivation above was performed for waves that propagate in vacuum, that is, within the light cone Q < uj/c 
and for real k. For evanescent waves with Q > uj/c and imaginary k the approach above has to be slightly modified, 
as it turns to be impossible to choose fictitious transmission amplitudes t'^'^ that guarantee energy conservation at 
the boundaries zi and Z2 of the fictitious cavity of Fig. [l] Nevertheless, the fictitious system may be easily altered to 
accommodate for evanescent waves [49] , and it turns out that the expression ( 24 ) and ( 25 ) remain valid even outside 
of the light cone [15]. Thus, the integration region of Eq. ([24]) may include real and imaginary values of k, as long as 
the wavenumber q = ui/c = \/ + fc^ is real, i.e., k should go along the imaginary axis from iQ to and then along 
the real axis towards infinity. 

The formalism developed above is a generalization to anisotropic systems of the formalism developed in Refs. |32| - 



[36] for the isotropic case. The resulting force (Eq. (24)) agrees with the usual Lifshitz's result expressed in terms of 
the reflection amplitudes in the isotropic case (r^^ — r^" = 0). 



The potential energy U of the system may be now obtained by integrating the force ( 24 ) with respect to the plate 
separation from oo towards the actual separation L, yielding 
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Notice that when L — > cx), A — 1 as e^*'"'^ — > for any positive of ry. We make a change of variable from k to q and 
perform the usual rotation in the complex plane from the positive real axis towards the imaginary axis to rewrite Eq. 
(1251) as 



du / d^Q log A 



U _ h 
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at r = 0, where we introduce an imaginary frequency w = iu, and as 
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(27) 



for r ^ 0, where we have accounted for the poles of f{iu) by performing residue- like integrations at the Matsubara 
frequencies ui — 2'KikBT/h and we define Ai = A{lu = iug). The prime in the summation means that the i — term 
should be divided by 2. Eq. ([26]) coincides with that derived in Refs. [30 ] [3T ] |42] for T = 0. 

Notice that for anisotropic plates, the matrices depend on their in-plane orientation. Therefore, the energy U 
(Eq. (25)) depends implicitly on the relative orientation 7 between the optical axes of the plates. Thus, we expect a 
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FIG. 2: Schematic diagram of the system, consisting of two parallel uniaxial non-magnetic plates whose optical axes lie on 
the surface. We indicate the optical axes of each plate with heavy double-headed arrows, and the projection of the axis of the 
second plate upon the first by a dashed double arrow. The z-axis is chosen to be orthogonal to the plates. We also indicate 
the plane of incidence. The angle ipi between the plane of incidence and the optical axis of plate 1 as well as the angle 7 from 
the optical axis of plate 1 towards that of plate 2 are indicated. The angle ^2 between the plane of incidence and the optical 
axis of plate 2 is ■!/)2 = — i^i — 7. The sign is due to the convention in Eq. (29 1. 



torque M on the plate s wh ich we may calculate simply by taking the derivative M — —dU/dj. It is easily verified 
that starting from Eq. (26 1 but setting Q = instead of performing the integral {A/An^) J dPQ yields Eq. (13) of Ref. 



[48j . i.e., the torque for a cavity in which the field is constrained to propagate along only one dimension, namely, along 
the normal to the surfaces. In ref. |48j the torque was obtained directly from the flux of angular momentum within 
the cavity. Here, we took a different approach, obtaining the torque from the angular dependence of the energy. The 
reason is that the angular momentum is well defined only for finite width beams and the overlap among the muUtiple 
reflections of a finite beam is incomplete for oblique incidence and ill defined for evanescent waves. 

We remark that our results are written in terms of the optical coefficient matrices Va of the walls of the cavity, 
about which we have made no assumptions. Thus, our results may be applied to systems of arbitrary absortance, 
conductivity and width, and they may be homogeneous or inhomogeneous. Up to this point the dependence of the 
energy, force and torque on the orientation of the plates has been implicit, through the unstated dependence of r^. 
In the next sections we will apply our result to specific cases where the dependence on 7 can be exhibited explicitly. 



III. SEMI-INFINITE UNIAXIAL SLABS 



We consider semi-infinite uniaxial non-magnetic crystals with their optical axes parallel to their surface. In this 
case, the reflection matrices ri and T2 may be obtained as particular instances for surfaces 1 and 2 of the formula [50] 

r = (82"^ +s^^cos6')"i(sj;^cos6l-S2"^), (28) 

which we derive in the appendix following the notation of Ref. 51 . Here, Si and S2 are the 2 x 2-matrices 
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where rio and rie are the ordinary and extraordinary refractive indices of the uniaxial crystal, 9 is the angle of incidence 
and ■0 is the angle from the plane of incidence to the optical axis of the crystal; its sign is chosen through the right 
hand rule around the normal of the surface that points outwards from the anisotropic medium. In Eqs. (29 1 we have 
assumed that the cavity is empty and thus we took its index of refraction of the cavity as 1. 

Consider the sketch depicted in Fig. [2] which displays the angle ipi of the optical axis of plate 1 with respect to the 
plane of incidence and the angle 7 between the optical axis of plate 2 and that of plate 1. According to our convention 
above, the angle -02 between the optical axis of plate 2 and the plane of incidence is -02 = — V'l — 7- We identify 
= V^ih "-o = V^j.! where ey and ej^ are the dielectric response functions of the plates along and perpendicular to 



their optical axes, respectively. The substitution of equations (29), (30) in (28) and (27 1 yields after a tedious algebra 
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FIG. 3: Casimir energy per unit area U /A between two ideal uniaxial plates as a function of the angle 7 between their optical 
axes for different temperatures T = 0,0.01,0.1,0.2,0.3,0.4,0.5 in units of Tic/LkB- 



an expression for the Casimir energy between anisotropic plates. We have checked that the resulting expression for 
A coincide with that obtained in Ref. [25] for the same system. Nevertheless, the calculation of Ref. [25] is directly 
applicable only to local uniaxial semi-infinite homogeneous systems, as their calculation is setup in terms of the bulk 
dielectric function of the plates. On the other hand, our result is written in terms of the reflection amplitudes 
of the plates, and thus may be applied to any system for which we can calculate these optical coefhcients, as shown 



explicitly in section IV 



Let us now consider an idealized case consisting of slabs which behave as perfect conductors along the optical axes 
but which are perfect insulators along the other principal directions. This could be realized through an array of 
perfectly conducting aligned wires insulated from each other. We further ignore the electric polarization across the 
wires, which would be appropriate, for example, if they were very thin. Then, each plate may be described by a 
dielectric response en — > 00 and ej^ 1. In this case the reflection matrix (28) becomes simply 



r(0,V): 



cos'^ 9 cos^ 'if) ~ sin 2-0 cos 0/2 



sin 2'0 cos6'/2 



— sin tj) 



which yields 



A = 1 ~ a^e^^''^ 



(31) 



(32) 



when substituted into Eq. (22), where 



= Tr(rir2) 



(cos 7 — sin^ 9 cos 0i cos ^2) 
(1 - sin^ ^cos^ 0i)(l - sin^ 0cos2 1/^2) ' 



(33) 



Eqs. (32) and (33) were previously obtained for this system in Ref. 
V'l + 7 



as becomes evident by writing the latter in 



terms oi ip = —ip2 

Substitution of Eq. ( 32 ) into Eq. ( 26 ) yields the T = energy of the system. Fig. |3] shows our numerical results 
for the energy of the idealized system discussed above. We have verified that they agree with the results presented in 
Ref. [27] where a path integral approach [55] [55] was used. Nevertheless, a similar substitution into Eq. (27) permits 



us to calculate also the T ^ energy of the system, shown in the same figure. Notice that for this idealized case 
there is a natural temperature scale Tic/kBL and a natural energy scale ATlc/L^, as they are respectively intensive 
and extensive quantities with respect to the area A of the interfaces. We use these scales to normalize the units in 
the figure. As T increases, the Casimir energy becomes more negative, so that the plates are more strongly bound. 
However, the slope of the energy as a function of the angle, and thus, the torque, diminishes with increasing T. This 
result may appear somewhat surprising, as it is well known that increasing the temperature yields an increase of the 
magnitude of the Casimir force. For high enough temperatures the energy becomes independent of the angle and 
approaches its value at 7 = 0, which is exactly half of the energy corresponding to t wo p erfect isotropic mirrors. This 
can be easily verified by substituting — 1 into Eq. (32) and the resulting A into (26) and (27). 



The torque for the ideal case may now be simply calculated by analytically deriving under the integral sign the 
energy (26) or (27) with respect to the relative angle 7 and performing the resulting integrals numerically. Fig. |4] 
shows that the torque M is a periodic function of 7 with period tt. It is null when the optical axes are aligned, 7 = 0, tt, 
corresponding to a stable equilibrium orientation. It is also null when they are orthogonal to each other, 7 = ±7r/2. 
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FIG. 4: Casimir torque per unit area M/A between two ideal uniaxial plates as a function of the angle 7 between their optical 
axes. The torque is calculated for the same temperatures as in Fig. Isl expressed in units of he/ Lks- 



corresponding to an unstable equilibrium. For T = the slope of M{^) seems singular at the stable equilibrium 
point. The existence of this singularity can actually be confirmed by deriving analytically the torque with respect to 
the angle 7 under the integral sign and examining the integrand after having taken the corresponding limit 7 — >■ 
and analytically performed the angular integral over ip. It turns out that the integrand has an infinite discontinuity 
at M = 0. This singularity had also been found previously in a ID calculation [48|, in which light is normally incident 
at the walls of the cavity. Furthermore, qualitatively the same behaviour of the torque as a function of the angle 7 
found in this ID calculation, is obtained in the 3D calculation for T = 0. The torque is not simply proportional to 
sin 27, so its extreme values are not at 7 = ±7r/4. However, as the temperature increases the torque becomes more 
sinusoidal-like. Note that the maximum torque per unit area is of the order of M/A 5 x 10"'^ x Tic/L^. Thus, for a 
separation L = 100 nm, the maximum torque per area unit is about 10~^N/m. It is interesting to compare this result 
with that of Rodrigues et al. [Jl] for the case of two slightly misaligned corrugated metals. It has been argued [41] 
that that system would yield the largest torques. Nevertheless, we obtained a torque of the same order of magnitude 
as their's. Furthermore, in our case the torque remains large over a wide angular range, independently of the area of 
the system, while in their case the torque is significant only within a very small angular range which decreases with 
the size of the system. 

Our theory may be applied to systems which are more realistic than the idealized case above. For example, in 
Fig. [5] we show the zero temperature torque between two uniaxial conductors whose response along the optical axis 
is characterized by the Drude dielectric function 



g||H = l- , ■ (34) 

For simplicity, we assume they do not respond along the perpendicular direction, ej^ = 1, as could correspond to 
parallel plates made up of an array of very thin aligned wires, each of which is a conductor with a finite Drude-like 
conductivity. We include results for different values of the electronic density and therefore of the plasma frequency 
Wp, as well as the ideal limit ojp — > 00. In the figure we used Wg = ttc/L as a natural frequency scale and we chose the 
electronic lifetime r = /uiq which would correspond to a typical values within a metal for L ~ lO^nm, although 
the results are very insensitive to r. In this case, wq would be of the order of a typical metallic plasma frequency. 

When ijjp 3> the plates behave as perfect mirrors and, as verified by Fig. [5] we recover the idealized case studied 
above, including the approach to the singularity in the slope at 7 = 0. As expected, the torque increases with lOp and 
becomes null at = 0. 



IV. FILMS 



Our formalism, summarized by Eqs. (26) and (27), is written entirely in terms of the reflection amplitudes ri and r2 
of the anisotropic plates through Eq. (22). This has the enormous advantage over previous formalisms [3T1[55J[2B] in 
that different systems may be explored simply by substituting their corresponding reflection amplitudes. For arbitrary 
uniaxial systems we could simply substitute the ordinary and extraordinary indices of refraction in the formulae derived 
in the previous section. In this section we further illustrate the versatility of our theory by calculating the torque 
between two anisotropic dielectric plates one of which is a film with finite thickness d. For the reflection amplitudes 
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FIG. 5: Casimir torque at T = between two uniaxial plates as a function of the angle 7 between their optical axes. The plates 
are conducting along their optical axis with a response ey described by a Drude dielectric function with a relaxation time r, 
and a plasma frequency LUp, and we assume the response perpendicular to the optical axis is e_L = 1. We chose a fixed value 
T = lOOO/oiQ, and we show results for various values of ujp = 1, 10 and 10^ x ujq where cjq = ttc/L. We also include the case 
ujp — 00 corresponding to the ideal case studied above. 



ri of the semi- infinite plate, which we take to be plate 1, we can simply use Eq. (28), but we have to calculate the 
reflection coefficients r2 of the film. 

Consider the reflection of an arbitrary polarized plane-wave with frequency lo from an uniaxial anisotropic free- 
standing homogeneous dielectric film of finite thickness d, with its optical axis parallel to its surface. We assume that 
the system is described by the same geometry as in the previous section, illustrated by Fig. [2j and that the film is 
non-magnetic. 

For each frequency lu and parallel projection Q of the wave-vector, two types of waves coexist inside the film, the 
ordinary and extraordinary waves, each of which has a wave vector component ±fc^ normal to the surface, where 
fi — o for ordinary waves, fi — e for extraordinary waves and we choose the upper sign for waves travelling or decaying 
towards z, i.e., we choose fc^ as the solution of the dispersion relation that is consistent with Im(kf^) > for evanescent 
waves and for propagating waves in the presence of finite or infinitesimal absorption. This gives a total of four waves 
in the film. Thus, we describe the electromagnetic field within the film as |51j 

F(z) = ^V^'^A'e''^'=''^ (35) 

a,iJ, 

where is a 4-component column vector which contains the electric and magnetic field projections parallel to the 
surface, F = {E^, By,Ey, —B^)^, 4"^ is the amplitude of the normal mode propagating in the az direction (cr = +, — ) 
with polarization /i = o, e, while is the eigenvector describing the electromagnetic field of a single mode with 
propagation direction a and polarization /i, and crfc^ ~ ±fc^ is the component of the corresponding wave-vector along 
z (see the appendix for details). 

The field reflected by the film has two contributions. One of them is the field reflected by its front surface. The 
other is the field transmitted back into the cavity from the film, after having entered the film from the cavity and 
being multiply reflected. Thus, we may write 

r2li2 = ro2^j2 + t2o'0- (36) 

where — {ipa^ '4'%)'^ is a- 2-component column vector containing the total o and e amplitudes of the waves travelling 
in the crz-direction within the film, = {^^2T^i2)^ the vector containing the p and s contributions of the incident 



field and is defined through Eqs. (17)-Jl8L the matrix Y2 is the sought reflection amplitude of the film, ro2 is the 



reflection amplitude of the front surface when light impinges from the cavity and therefore it is given by Eq. ( 28 1 , and 
t2o is the transmission amplitude from the film into the cavity across the front surface. Similarly, the field that travels 
towards z within the film is the sum of the field transmitted from the cavity into the film plus the field travelling 
towards —z and internally reflected back at the front surface, so we may write 

^+ = to2^j2 + r2o'0-, (37) 



where to2 is the transmission amplitude from the cavity into the film through the front surface, and r2o is the internal 
reflection amplitude of the front surface when light impinges from within the film. Similarly, at the rear surface, we 
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FIG. 6: Casimir torque at T = between a semi-infinite uniaxial, dissipation-less conducting slab and a film made up of the 
same material as a function of thickness d of the film. The slabs are modelled as Drude metals with a plasma frequency Wpy 
along the optical axis twice as large than its value along the perpendicular direction. The curves are labelled by the value 
of Up^/ujo — 1/4, 1/3, 1/2, 7/10, 1, 2, 3, 5, where uuo = nc/L. We fixed the relative angle between the optical axes at 7 = 7r/4. 



may write 

K~^xjj^ = V2I3K,'4>+ (38) 

where k = diag(e''^°'^, e''^^'') is a 2 x 2 diagonal matrix that accounts for the phase acquired by the ordinary and 
extraordinary waves as they travel from the front surface towards the rear surface of the film. Notice that the internal 
reflection matrix of the rear interface coincides with the internal reflection matrix at the front interface as we consider 



a free standing fllm, and hence we denote both matrices by r2o- Elimination of ■0+ and ■0_ from Eqs. (36l-(38) yields 



r2 = i"02 + t2o/tr2o/«(I - (r2o'«)^) ^to2. (39) 

We remark that the same result would be obtained by generalizing Airy's method [52], summing the amplitudes of 
successive multiple reflections and refractions, but taking account of their 2 x 2-tensorial character due to the mixing 
of polarizations at the interfaces. Equivalent results may also be obtained through a 4 x 4 transfer matrix formalism 
[51LI53) . though care must be taken to avoid numerical instabilities. The calculation of the required matrices to2, t2o, 
and r2o, as well as the calculation of the wave vector components is shown in the Appendix. 



By substituting Eq. (39 1 into Eqs. (26 1 or (27) we finally obtain the Casimir energy corresponding to an anisotropic 
film in front of a semi-infinite, anisotropic substrate. Then, the torque is obtained as before, simply by deriving 
analytically the energy under the integral sign with respect to the angle 7 between the optical axes. 

In Fig. [6] we show the torque between a semi- infinite anisotropic substrate and an anisotropic fllm as a function of 
the thickness d of the fllm. The substrate and the film are made of the same material which we assume is an uniaxial 
conductor with different plasma frequencies along its principal directions, i.e., with a tensorial anisotropic effective 
mass, and which we model by a Drude-like response along and perpendicular to the optical axis 

ep(^) = l- 2!'" / ' A^=^'ll- (40) 

LU^ + lUJ/T 

We include results for different values of ujp± for a fixed quotient iOp^\/ujp± ~ 2 and for simplicity, we disregarded the 
dissipation, assuming r = 00. 

From Fig. [6] we see that for very thin films the torque becomes smaller as the normalized film thickness d becomes 
smaller and eventually the torque becomes null for a zero film thickness d — 0, as r2 — vanishes in that limit, as 



can be verified from Eq. ( 39 1 . For large enough values of d the torque tends to its asymptotic value corresponding to 
semi-infinite plates. Curiously, the figure shows optimal values of the film thickness for which the torque is maximized. 
For the chosen parameters, the maximum torque may be as large as twice its asymptotic value. The optimal thickness 
shifts towards smaller values as the plasma frequency ujp± is increased. Furthermore, the magnitude of the torque at 
the optimal thickness is maximum for values of ujpj_ around ojq. 

In order to explain qualitatively the results shown in Fig. [6| we consider a ID model in which light is normally 
incident at the walls of the cavity. In this case, it is possible to decouple the ordinary and extraordinary rays inside 
the media by choosing the polarization of the incident electrical field. Light polarized in the directions normal or 
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FIG. 7: Magnitude of the normal-incidence reflection coefficients [rj | of an anisotropic, thin, dissipation-less Drude conducting 
film with plasma frequencies cjpy = 2ujp^, and magnitude of the reflection anisotropy |Ar2| as a function of the film thickness 
d. The frequency is lj = Q.lu}y±. 



parallel to the optical axis will only excite the ordinary or extraordinary wave inside the media and will be reflected 
without changing its polarization with reflection amplitudes r"*" or respectively. The torque at T = for this one 
dimensional system can be calculated using the results of Ref. [IS] , 

;isin27 [°° AriAr2e~2«V' 



27r 



M = / du ^ , (41) 



AriAr2 sin^ ^e^^uL/c .)_ (i _ ^\Vf.\}^^-2uL j c-^{^^ _ ^j_^j_g-2«L/c-) 



where Ar^ = r a — r-^ . The reflection coefficients are 



(42) 



for plate 1, and 



e 



2ik,,d 



'2 " 1 - (r5')2e2 



(43) 



for the film, where /i =||, -L and = Jt^Lojc is the wave number inside the film corresponding to /^-polarization. 



According to Eq. (41 1, we might understand the dependence of the torque M on the thickness d of the film if 
we first understand the dependence of the anisotropy Ar2 on d. Consider first a wave of frequency to > a;p_L,a;p||, 
normally incident on the film 2. In this case, the film would be transparent for both polarizations and the anisotropy 
would be small. For intermediate frequencies, ujp± < cj < Wp|| , the film would be a good reflector for one polarization 
and transparent for the perpendicular polarization, yielding a large anisotropy which nevertheless would be quite 
insensitive to the thickness of the film. However, for uj < ijjp± < Wpy the penetration depth £± for polarization 
perpendicular to the optical axis would be larger than the penetration depth for polarization along the optical 
axis. For a sufficiently thick film, d > £± > and the film would be a good reflector for both polarizations as no 
radiation would get across it, so the anisotropy would be small. For a sufficiently thin film i± > > d, the radiation 
would reach the back surface and leave the film, which would be a poor reflector for both polarizations and thus 
the anisotropy would again be small. In the limit of vanishing thickness d — > all the incident radiation would be 
transmitted through the film and none reflected. Nevertheless, for an intermediate thickness such that £± > d > 
the film would be a good reflector for polarization along the optical axis but a poor reflector along the perpendicular 
direction, yielding a large anisotropy. 

These qualitative features are displayed in Fig. [Tj where we show the magnitude of the reflection coefficients {rl^l as 
well as the magnitude of the anisotropy |Ar2| as a function of d for the case uj = 0.1ujp±. In the limit of large thickness 
d S> c/ujp±, the anisotropy of the film Ar2 tends to its asymptotic value, which is non zero due to the relative phase 
between the complex optical coefficients corresponding to both polarizations. The film thickness which maximizes 



the anisotropy can be calculated from Eq. ( 43 1 in a straightforward way, yielding 



|c^l<^P±,Wpy (44) 

Ldp±UJp\\ 
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FIG. 8: Casimir torque ID (only light normally incident is accounted for) at T = of the same system considered in Fig. [b] 
as a function of the thickness of the film d/L. The torque is calculated for cjpx = ujp = 1/10, 1/7, 1/3, 1/2, 1, 2, 3, 5 (in units of 
cjo = nc/ L). The rest of the parameters are the same as in Fig. |6] 
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FIG. 9: Casimir torque at T = between a semi-infinite substrate of TiOs and a film of calcite, as a function of the thickness 
of the film d. The torque is calculated for several separations L — 10, 50, 100, 200, 400, 600, 800, 1000 nm. The angle between 
the optical axes is 7 = 7r/4. 



for small real or imaginary frequency. 

In Fig. [8] we show the torque as a function of the width d for the same system as in Fig. [6] but assuming that 
the field propagates only in ID [15]. The results display the same qualitative behavior as in the 3D case, although 
the plasma frequency ujp± for which the magnitude of the torque attains its maximum is shifted to a smaller value 
(jjpj_/ujQ « 1/3 compared to 7/10 in Fig. [6j We may estimate the optimal width for each value of ujpj_ through 
the following considerations: The exponential in Eq. (41) suppresses the contribution of large imaginary frequencies 
u > c/2L to the Casimir effect. Thus, in the retarded regime LOp± > coq we may consider mostly low imaginary 
frequencies u/ujp^ < 1 to understand the effect. Furthermore, as the relevant frequency scale is c/L ^ ujq/tt, then the 
localization of the optimal film thickness in the retarded regime could be estimated by substituting |a;| = u = ujo/tt 
within Eq. (44). Thus, we obtain d™/L « 0.025, 0.01 and 0.004 for the choices u!p±/ijjQ = 2,3 and 5 respectively. 
These estimates are in good agreement with the ID calculation. Fig. [8j In the 3D case one would have to take into 
account modes that propagate along non-normal directions. However, as the reflection amplitude increases towards 
unity for any polarization as the angle of incidence moves away from the normal direction, the anisotropy in the 
reflectance and the corresponding contribution to the torque also diminish. Thus, the most important contributions 
to the torque come from modes that propagate close to the normal and our estimates for are also in good qualitative 
agreement with our full 3D calculation. Fig. |6] 

The simple estimates above fail in the non-retarded regime ujp± < ujq for both the ID and 3D calculation, as in 
this case a wider range of frequencies, going beyond a;p_L,a;p||, contributes to the torque. Furthermore, for small 
separations L the Casimir effect in the 3D case is dominated by surface plasmons [SS] which are completely left out 
of the ID calculation and of our estimates above. 

To conclude this section, in Fig. [9]we show the Casimir torque between a crystal of BaTiOa and a crystal of calcite. 
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as in Ref. [53]. However, instead of two semi-infinite crystals, we consider a thin calcite film over a semi-infinite 
BaTiOs substrate with vacuum in between and we vary the thickness d of the film. We use the same model as in Ref. 
[23] to describe the dielectric response of the plates along each of their principal directions, namely, two undamped 
oscillators [55] to account for both the IR and UV resonances. We show results for T = which we expect to hold 
even at room temperature for separations in the range L < 1/im. We have confirmed that our numerical results, 
performed at T = 0, are in agreement with those of Ref. [M] in the limit of very thick plates. 

As expected. Fig. [9] shows that for very thin films the torque becomes null, while for thick enough films it attains 
its asymptotic value corresponding to semi-infinite plates. Fig. [9] displays small maxima corresponding to an optimal 
thickness for which the torque is maximized for each value of the separation, analogous to those found above for 
anisotropic conducting films, although these maxima are very broad. In this case, the maximum torque is about one 
order of magnitude smaller than for the case of anisotropic conductors shown in Fig. [6] when typical metallic plasma 
frequencies are chosen, and about two orders of magnitude smaller than for the case of the ideal uniaxial plates shown 
in Fig. [4] 

In order to observe the Casimir torque, several experimental setups have been suggested (Ml [57]. If the vacuum 
Casimir cavity is replaced by a cavity filled with a dielectric liquid, and if the dielectric function of the fluid is 
intermediate between that of the cavity walls, then the sign of the Casimir force may be reversed, becoming repulsive 
instead of attractive. Munday et al. [24] took advantage of this sign reversal and estimated that a quartz or calcite 
anisotropic disk with a diameter and a width of some tens of fim would float at a height of about one hundred nm 
over the surface of a barium titanate anisotropic crystal if immersed in ethanol. At this height, the repulsive Casimir 
force would balance the weight of the disk, which would then be able to rotate freely around its axis. Then it would 
be possible to first align the disk with the polarization of a laser beam and optically monitor the rotation towards 
the equilibrium orientation after the laser beam is turned off. The characteriztic time would depend of the viscosity 
of the fiuid and on the magnitude of the driving Casimir torque which would thus be obtained. 

In a second proposal by the same group [ST] a large set of much smaller disks of /zm scale width and diameter are 
kept separated a few nanometers from the substrate not by a repulsive Casimir force but by surfactant molecules that 
attach to the surfaces of both the substrate and the disks which in this case would be immersed in an electrolyte. 
The orientation of such small disks would be disrupted by their rotational Brownian motion, but the distribution of 
their orientations would yield **** Sepasuchi. 

According to our calculations above, structured materials made up of aligned nanowires embedded in dielectric 
matrices are subject to much larger torques than the birefringent crystals considered my Munday et al. [241 157 ] . 
Thus, the Casimir torque could be measured using a torsion balance. To avoid alignment problems, a spherical 
anisotropic particle could be used instead of a fiat disk. Hanging a 50 /im sphere from a 10 fim long aluminum wire 
of 50 nm radius 

frecuencia =100Hz vacfo placa giratoria. 50HZ. 

V. CONCLUSIONS 

We have extended a fictitious-cavity approach to calculate the Casimir force, energy and torque between anisotropic 
media with a planar geometry in terms only of their optical coefficients. Our results are applicable to arbitrary 
anisotropic materials which may be conducting or insulating, opaque or transparent, semi-infinite or with a finite 
thickness, with a non-dispersive or a frequency dependent response, homogeneous or structured, as long as we can cal- 
culate their corresponding reflection coefficients. Our expressions for the torque were simply obtained by analytically 
deriving the energy with respect to the relative orientation of the plates. The formalism has allowed us to perform 
calculations at zero and at finite temperature. 

We have reproduced some known results and we have calculated the torque for ideal uniaxial systems which are 
the anisotropic counterparts to the ideal Casimir mirrors, namely, systems which are perfect conductors along some 
directions and perfect insulators along others. We also calculated the torque for systems with a Drude conductivity. 
As a non trivial application of our formalism, we have calculated the torque for a system consisting of two anisotropic 
slabs, one of which is a semi-infinite substrate and the other a thin film of thickness d. The slabs were modelled as 
Drude metals and made up of the same material in order to study the effect of the film thickness in the torque. Our 
numerical results show that exists an optimal value of the film thickness where the torque is maximized. We have 
estimated the optimal thickness in the retarded regime by studying the anisotropy of the skin depth evaluated at 
the characteristic frequency of the cavity; the maximum torque is achieved for films whose thickness is intermediate 
between the largest and the smallest of the two principal skin depths. Finally, we calculated the Casimir torque for 
a thin film of calcite above a barium titanate substrate. The results seem qualitatively similar although the torque is 
smaller than that for uniaxial conductors and for the ideal system. 

In conclusion, we developed a very general formalism to calculate the Casimir torque and illustrated its use by 
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studying the torque at zero and finite temperatures for ideal and realistic, conducting and insulating systems of semi- 
infinite and finite widths. We expect that the simplicity and generality of our formalism motivates further studies of 
these relatively unexplored aspects of the Casimir effect between anisotropic media. 
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Appendix 



In this appendix we derive the expressions for the matrices t2o, r2o, to2 and the ^-component of the wave- vectors 
of the ordinary and extraordinary waves. We derive these expressions by following the method developed in Refs. 
[51] [53]. We assume that the system is described by the geometry shown in Fig. |2] Thus, the dielectric tensor is 
given by 



ey cos^ ■(/' -|- e_L sin^ V (e_L — ey) sinV'COS V' 
(ej- — ey) sin?/'cos?/' ej_ cos^ V + ey sin^ V' 
e_L 



(45) 



All of the quantities above (the dielectric functions and the angles) refer to plate 2, but we omit the corresponding 
index (2) in order to simplify our notation. On the other hand, from Maxwell equations, the components of the 
electromagnetic field parallel to the x — y plane, can be cast in a set of four differential equations 



where 



d,F = iqTF, 



F — {Ex, By, Ey, —Bx)^, 



(46) 



(47) 



Ex, Ey, Bx and By are the electric and magnetic field components parallel to the interfaces, q is the free-space 
wavenumber and F is a 4 x 4 matrix with elements 



rii 


= -sm0e^x/<^zz, 


ri2 


= 1 - sin^6'/ezj;. 


ri3 


= -sm9e^y/e^^, 


ri4 


~ ^24 = ^31 = ^32 


r2i 


— ^xx ^xz^zx / ^zz , 


r22 


= - sin Ocxz/ezz, 


r23 


— ^xy ~ ^xz^zy / ^zzj 


r34 


= 1, 


r4i 


— ^yx ^yz^zx 1 ^zz-j 


r42 




r43 


= tyy - Sin^ 9 - tyz 



33 



(48) 



^zy/^zz- 



where 9 is the angle of incidence. Within a homogeneous medium T is independent of z and Eq. (46) has four 
particular solutions of the form 



ikpz 



1,2,3,4, 



(49) 



where kg is the component of the propagation vector parallel to the z-axis. Substitution of Eq. (491 into the Eq. (46) 
yields the eigenvalue equation 



(M-gr)F,(o) = 



(50) 



whose eigenvalues 



(51) 
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are 



ko = Jq, 



(52) 



where the quantities / and J are defined by Eqs. ( |30[ ). The corresponding four eigenvectors are 

■ cot jj \ 

VI 



(53) 



± tan?/)/ J 
ej_ tan -0 / 
±1/J 
1 



— cot ll) 

1 



(54) 



where the eigenvector V^ describes the /i-polarized wave with a wave vector component ak^ aiong the z-axis. 

Consider the refiection of light impinging at the interface between vacuum and a semi-infinite anisotropic medium 
from within the medium. The electromagnetic field at the vacuum side is made up of only a transmitted wave. Inside 
the medium, the electromagnetic field is made up of four waves F — Xicr /i V'^K^s*''''^^- The field F is continuous 
across the interface, thus 



Ft = V V^V;^e*'"'=^'^ 



(55) 



Writing the generalized vector field Ft in terms of the spinor ^t = (Cf,Ct)'^, defined through Eqs. (n7|-(n8|, when 
they are applied to the transmitted field, we get 



Ft ^i~^Pcos9,^let-^t cos 



(56) 



Substitution of this equation into Eq. ( 55 ) yields four linear algebraic equations that can be recast as two matrix 
equations for the spinors ^t and V± = (V'±'V'±)^j 



and 



where Si and S2 are the 2 x 2-matrices 



Sl 



cos 9^t = 81-0+ + S3l/>_ 



»2 I -sro -tre 



(57) 
(58) 

(59) 



here, V/^^ is the i-th component of the eigenvector given by Eqs. (54) and when these are substituted in Eqs. (59) 



yield the same expressions than those of Eqs. (29). Similarly, S3 and S4 are the matrices 



S3 = 



S4 = 



yi- yl- 



Elimination of from Eqs. (57)-(58l yields 

■0+ = — (si + S2 cos &)^^{&-^ + S4 cos G)"^- 



Thus, from = V2q'4>- we obtain 



I"20 = -(Sl + S2 COS0) ^(S3+S4C0S( 



The matrix can be obtained by substituting Eq. (61 ) into Eq. (58) to get 

it = (S4 + S2r2o)V'-, 

from which we obtain 

t20 = (S4 + S2r2o). 



(60) 

(61) 
(62) 
(63) 
(64) 
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Finally, in order to calculate the matrices to2 and ro2 , we consider the reflection of light impinging at a semi-infinite 
anisotropic medium from vacuum. The electromagnetic field at the medium side is made up of two transmitted waves, 
F — J2fj, V'+^^e*'^f ^. Imposing continuity on F across the interface 



(65) 



Writing the generahzed vector field Fi and Fr in terms of the spinors — (^f,^f)^ and — {£,r,£,r)'^ we get 

F,^(ercos0,ef,ef,e^os0)^ (66) 

and 

Fr = {~er COS e, COS Of. (67) 



Substitution of this equation into Eq. ( 65 ) yields four linear algebraic equations that can be recast as two matrix 
equations for the spinors and 



cos0{i, - ir) = Sl'/'^ 



and 



Elimination of from the last two equations yields 



where ro2 is given by Eq. (28 1. Elimination of from Eqs. (68)-(69) yields 

to2 =Sj;^cos0(I-ro2). 



(68) 
(69) 
(70) 
(71) 



This concludes the calculation of the reflection and transmission matrices corresponding to the front surface of the 
film. 
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